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Abstract 

A method for describing the quantum kink states in the semi-classical limit of several 
(l-l-l)-dimensional field theoretical models is developed. We use the generalized zeta function 
regularization method to compute the one-loop quantum correction to the masses of the 
kink in the sine-Gordon and cubic sinh-Gordon models and another two P{(j))2 systems with 
polynomial self-interactions. 
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1 Introduction 



BPS states arising both in extended supersymmetric gauge theories ^ and string/M theory 
play a crucial role in the understanding of the dualities between the different regimes of the system. 
In this framework, domain walls appear as extended states in A/" = 1 SUSY gluodynamics and 
the Wess-Zumino model 0. This circumstance prompted the question of whether or not these 
topological defects saturate the quantum Bogomolny bound. A return to the study of quantum 
corrections to the masses of (l+l)-dimensional solitons has thus been unavoidable. These subtle 
matters were first addressed in the classical papers of Dashen, Hasslacher and Neveu, and 
Faddeev and Korepin, 0, for the purely bosonic \[(p]2 and sine-Gordon theories, and then in 
Reference for the super-symmetric extension of these theories. Analysis of the ultraviolet cut- 
off regularization procedure in the presence of a background is the main concern in the papers of 
Reference [0 : the authors carefully distinguish between using a cut-off either in the energy or in the 
number of modes. The second method leads to the same result as in the computation performed 
by DHN for bosonic fluctuations. Another point of view is taken in Reference [|], where SUSY 
boundary conditions related more to infrared behaviour, are carefully chosen. On this basis, and 
by using higher-derivative ultraviolet regularization (SUSY preserving), the authors demonstrated 
an anomaly in the central charge that compensates for the extra (quantum) contribution to the 
classical mass. 

In this paper we shall confine ourselves to purely bosonic theories and leave the treatment of 
fermionic fluctuations for future research. We address the quantization of non-linear waves relying 
on the generalized zeta function regularization method to control the infinite quantities arising 
in the quantum theory. This procedure has been used previously in computing Casimir energies 



and the quantum corrections to kink masses, see ||TT[. We shall develop this topic, however, in a 
completely general way, also offering a comparison with other approaches. As well as obtaining 
exact results, we also shall explain how asymptotic methods lead to a very good approximation 
of the right answer. We believe that the novel application of the asymptotic method should be 
very useful in the cubic sinh-Gordon model as well as in multi-component scalar field theory, 
where the traditional approach is limited by the lack of detailed knowledge of the spectrum of the 
second-order fluctuation operator (see []13| , [|T^ for extensive work on multi-component kinks and 
their stability) . 

The organization of the paper is as follows: In Section §.2 the general semi-classical picture 
of quantum solitons, the zeta function regularization procedure, the zero-point energy and mass 
renormalization prescriptions, and the asymptotic method are described. In Section §.3, we apply 
the method to the "loop" kinks of the sine-Gordon , A(0)2, and cubic sinh-Gordon models. In 
the first two paradigmatic cases, it is possible to obtain an exact result, which allows comparison 
with other methods. Approximate computations by means of the asymptotic expansion of the 
heat function are also offered to test the goodness of our procedure against the well known exact 
answers. Section §.4 is devoted to the analysis of the "link" kink arising in the A(0)2 model. 
Finally, Section §.5 offers an outlook on further applications of our approach. 
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2 Semi-classical picture of quantum soliton states 



We shall consider (H-l)-dimensional scalar field theories whose classical dynamics is governed by 
the action 



1 dip dtp 



2 dyi'dyf. 

We choose the metric tensor in T^(R^'^) as g = diag (1, —1) and the Einstein convention will be 
used throughout the paper. We shall not use a natural unit system because we wish to keep 
track of h in our formulas; nevertheless, we choose the speed of light to be c = 1. Elementary 
dimensional analysis tells us that [h] = ML, [[/{ip)] = ML~^ and [ip] = M2L2 are the dimensions 
of the important quantities. 

The classical configuration space C is formed by the static configurations ip{y), for which the 
energy functional 



is finite: C = {ip{y) / E{ip) < +00}. In the Schrodinger picture, quantum evolution is ruled by 
the functional differential Schrodinger equation: 

th^^^^{y),t]=Il^^{y),t]. 



The quantum Hamiltonian operator 



dy 



5 



2 Sijiy) 5i,{y) 



+ E 



acts on wave functionals '^[^{y),t\ that belong to L'^{C). 

We wish to compute the matrix element of the evolution operator in the "field" representation 



for the choice 



V[^{y,t)]exp{--S[^] 



i;^'\y,o) = My) , i^^f\y,T) = My) 

where ipK{y) is a kink static solution of the classical field equations. We are, however, only 
interested in the loop (h) expansion of G up to the first quantum correction. Also performing a 
analytic continuation to "Euclidean" time, t = —ir, T = —ipi, this is achieved by the steepest- 
descent method applied to the Feynman integral in (|1|): 



Ge {ipK{y), ^k(i/); 13) = exp 
where /C is the differential operator 



h 



Det"2 



P/C 



:i+oih)) 



+ 



dy^ dip"^ 



3 



and P is the projector over the strictly positive part of the spectrum of KL. Note that, on the math- 
ematical side, the steepest-descent method provides a well defined approximation to the Feynman 
integral if the spectrum of the quadratic form /C is positive definite and, on the physical side, 
the zero eigenvalue that appears in Spec(/C) contributes to the next order in the loop expansion: 
it is due to neutral equilibrium on the orbit of the kink solution under the action of the spatial 
translation group. Moreover, in order to avoid the problems that arise in connection with the 
existence of a continuous spectrum of /C, we place the system in a interval of finite but very large 
length L, i.e. x e [— f , f ], and, eventually -after assuming periodic boundary conditions on the 
small fluctuations all throughout the paper- we shall let L go to infinity. 
^From the spectral resolution of /C, 

l^Uv) = ^lUv) , ojI e Spec(/C) = Spec(P/C) + {0}, 
we write the functional determinant in the form 

^2 



Det 



+ IC 



JJdet 



Q^2 + 



All the determinants in the infinite product correspond to harmonic oscillators of frequency io;^ 
and thus, with an appropriate normalization, we obtain for large (3 



^'n(S) 



■E„^n(l+o(ft)) 



where the eigenvalue in the kernel of IC has been excluded. 
Inserting eigen-energy wave functionals 

we have an alternative expression for Ge for /3 — > oo: 

Ge {My),Myy,P) = %[My)] ^^oiMy)] e 

and, therefore, we obtain 



a;2>0 

PIC 



|*o[V'k(?/)]| =Deti 



TT 



+ o{n') , 



as the kink ground state energy and wave functional up to one-loop order. 
We define the generalized zeta function 



Cpic(s)=IV(P/C)- 

associated to the differential operator P/C. Then, 
4 = -BlV'Kl + 5Tr(PK)3+o(/i2) 



1 



UJ. 



2\s 



(2) 
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\MMy)]f = vrnexp{-Cpyc(0)}exp |--A^(0)| (3) 

show how to read the energy and wave functional of the quantum kink ground state in terms of the 
generahzed zeta function of the projection of the second variation operator /C in the semi-classical 
limit. 



2.1 The generalized zeta function regularization method : zero-point 
energy and mass renormahzations 

The eigen-functions of JC form a basis for the quantum fluctuations around the kink background. 
Therefore, the sum of the associated zero-point energies encoded in Cpk:(— |) in formula (||) is 
infinite and we need to use some regularization procedure. We shall regularize Cp/c(— |) by defining 
the analogous quantity Cpic{s) at some point in the s complex plane where Cpic{s) does not have a 
pole. Cpic{s) is a meromorphic function of s, such that its residues and poles can be derived from 
heat kernel methods, see [1^. If Kic{y, z; (3) is the kernel of the heat equation associated with /C, 

^ +K\K^{y,z-(3) = Q , K^{y,z-Q)=5{y-z) (4) 



dl3 

the Mellin transformation tells us that. 







where, 

hp^[[5] = Tre-^^'^ = Tre'^'^ - 1 = -1 + j dy Kic{y,y; P) 

is the heat function hpjc[(3], if /C is positive semi-definite and dimKer(/C) = 1. Thus, the "regu- 
larized" kink energy is in the semi-classical limit: 

e^{s) = E[<j>K] + ^f^''^\pKis) + oih') (5) 

where is a unit of length"^ dimension, introduced to make the terms in @ homogeneous from 
a dimensional point of view. The infiniteness of the bare quantum energy is seen here in the pole 
that the zeta function develops for s = — |. 

To renormalize we must: A. Subtract the regularized vacuum quantum energy. B. Add 
counter-terms that will modify the bare masses of the fundamental quanta, also regularized by 
means of the generalized zeta function. C. Take the appropriate limits. 

A. The quantum fluctuations around the vacuum are governed by the Schrodinger operator: 



rf2 d^U 



dy"^ d^p'^^ 



where V'v is a constant minimum of ?7[V']- The kernel of the heat equation 

-^ + v]Kviy,z;P) = , Kviy,z;0) = Siy-z) 



5 



provides the heat function hv{P), 



/iv(/9)=Tre-^^ = j dy Kv{y,y; P). 
We exclude the constant mode and, through the Melhn transformation, we obtain 

Cv{s) = -— / dpp'-'Tie-^''. 

r(s) Jo 

The regularized kink energy measured with respect to the regularized vacuum energy is thus: 

e^^is) = E[0k] + Ai£K(s) + o(;i2) 

= E[M + ^f^''^' [(pds) - (vis)] + o{h'). 



B. If we now go to the physical limit = lim^^_ 



we still obtain an infinite result. 



The reason for this is that the physical parameters of the theory have not been renormalized. 
It is well known that in (l+l)-dimensional scalar field theory normal ordering takes care of all 
the infinities in the system: the only ultraviolet divergences that occur in perturbation theory 
come from graphs that contain a closed loop consisting of a single internal line, [|16] . From Wick's 
theorem, adapted to contractions of two fields at the same point in space-time, we see that normal 
ordering adds the mass renormalization counter-term. 



dy 5m^ 



d^U 



: +o{h^ 



to the Hamiltonian up to one-loop order. To regularize 

^ . dk 1 

dm 



47r ^k^ + U"{^y) 

we first place the system in a ID box of length L so that 6m^ = j^Cvij), if the constant eigen- 
function of V is not included in (^y. Then, we again use the zeta function regularization method 
and define: Sm'^{s) = ~x^^^j^Cv{s + 1). Note that drn^ = lim^^_i 5m?{s). The criterion behind 
this regularization prescription is the vanishing tadpole condition, which is shown in Appendix B 
of Reference to be equivalent to the heat kernel subtraction scheme. 
The one-loop correction to the kink energy due to H{5m?{s)) is thus 



K/ 



(tAkI H{5m\s)) IV^k) - (^vl H{5m\s)) |^v) 



2s+l 



lim — /i 

L— >oo 2ij 

— lim 

L— >00 2i7 



T{s) 
Tis + 



dyCv{s + 1) 



'd^U 




d^U 




d^^ 


V'K 




V'v. 



2s+l 



Tis) 



-Cvis + 1) / dyViy) 



(6) 



because the expectation values of normal ordered operators in coherent states are the correspond- 
ing c-number-valued functions. 
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C. The renormalized kink energy is thus 



(7) 



whereas the renormahzed wave functional reads 

PIC 



Det4 



Det~ 



V 



TThexp {- (Cpyc(O) - Cv(0))} exp <; -- 



<^(0) - ^(0) 
as as 



2.2 Asymptotic approximation to semi-classical kink masses 

In order to use the asymptotic expansion of the generahzed zeta function of the JC operator to 
compute the semi-classical expansion of the corresponding quantum kink mass, it is convenient 
to use non-dimensional variables. We define non-dimensional space-time coordinates = rridy^ 
and field amplitudes (f){x^) = Cdipiy'^), where and Cd are constants with dimensions [m^] = 

11 — 

and [cd] = M~2L~^ to be determined in each specific model. Also, we write Uicp) = -^Uiip). 
The action and the energy can now be written in terms of their non-dimensional counterparts: 



Em 



c 



The important point is that the Hessians at the vacuum and kink configurations now read 



V = mi 



K, = m„ 



dx' 



where 



and 



— V[x). Therefore, 



m 



2s 



Cvis) 



CM 



The asymptotic expansion is superfluous if Tre"'^^'^ and (plci^) are susceptible of an exact 
computation. If V{x) is a potential well of the Posch- Teller type, see [1^, one can completely 
solve the spectral problem for K. and there is no need for any approximation to Cpici^)- general 
the spectrum of K, is not known in full detail, specially in systems with multi-component kinks, 
and we can only determine Cpjc{s) by means of an asymptotic expansion. Nevertheless, we shall 
also compute the asymptotic expansion of Cpjci^) ^^e cases where the exact answer is known in 
order to estimate the error accepted in this approach. 

In the formulas (^), (|^), (||) and (0) we replace V, /C and v"^ by V, IC and v"^ and write the 
kernel of the heat equation for K, in the form: 



K^{x,x'-(3) = K^{x,x'-(3)A{x,x'-(3) 
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A{x,x'; (3) is thus the solution of the PDE 

X — x' d d'^ 



d_ 

d(3 



(3 dx dx^ 



V{x) A{x,x';(3) = 



with "initial" condition: A{x,x';0) = 1. 

For /5 < 1, we solve (||) by means of an asymptotic (high-temperature) expansion: A{x, x'\ /3) = 
^"(•^' -^O/^"- Note that there are no half-integer powers of (5 in this expansion because our 
choice of boundary conditions with no boundary effects. 

In this regime the heat function is given by: 



Tr e 



2 



dx Kj^{x, x; P) 



;2 OO 



n=0 



/ dxan{x,x)(3" 

2 



n=0 



It is not difficult to find the coefficients a„(x, x) by an iterative procedure starting from ao(a;, x') = 
1. This procedure is explained in the Appendix, which also includes the explicit expression of some 
of the lower-order coefficients. 

The use of the power expansion of hp^(5\ = TTe~^^^ in the formula for the quantum kink 
mass is quite involved: 

1. First, we write the generalized zeta function of V in the form: 



Cvis) 



1 mgL 
T(s) 



d/3(3'—^e 



+ By(s) 



with 



Bv(s) 



mdLT[s - 

4^ v^'-^r[s] 



Cvis) 



47r r s) 



and r[s,f^] and 7[s — being respectively the upper and lower incomplete gamma functions. 



see 



r8| . It follows that Cv{s) is a meromorphic function of s with poles at the poles of 7[s — |, f 



which occur when s — | is a negative integer or zero. B^{s), however, is a entire function of s. 

2. Second, from the asymptotic expansion of hj^[j3] we estimate the generalized zeta function 
of PJC: 



(pds) 



T{s) 



■'^ V 47r ^^^^ Jo 



+ 



sT{s) Tis)^ 



2^ an(/C)- 



2{s+n~\) 



n<no 



where 



2^ a„(/C)- 



n>no 



2{s+n-i) 



is holomorphic for Re s > —uq + ^, whereas 



dl3Tie-^^^p'-^ 
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is a entire function of s. The values of s where s + n — | is a negative integer or zero are the poles 
of Cpk.{^) because the poles of 7(5 + n — 1,^^] lie at these points in the s-complex plane. 

Renormalization of the zero-point energy requires the subtraction of Cvi^) from Cpici^)- We 
find, 



1 



Vis) 



no-l 



n=l 



47r 



2(s+n-i) 



and the error in this approximation with respect to the exact result to AiS^ is: 

hmd, 



Note that the subtraction of Cvl-s) exactly cancels the contribution of ao{)C) and hence, the diver- 
gence arising ect s — n — 0. The quadratic ultraviolet divergences appear in this scheme as 
related to the pole of Cvi^) at s = — |, n = 0. 

3. Third, Ai£^ now reads 



Ai£ 



K 



hnid , h 

+ — hm 



2s 

no-l 

E 

n=2 



ml 



47rr(s) V 



2s+l 



4vrr(-i) 



^2n-2 



The logarithmic ultraviolet divergences, hidden at first sight in the DHN approach, arise here in 
connection with the pole of Cp)c{^) ~ Cv('^) Bi, s — n — 1. 

The surplus in energy due to the mass renormalization counter-term is. 



hai{IC) 
lim — hm 

= -z^aiiJC) lim 

and the deviation from the exact result is 



\mdj T{s) 



rrid 



2s+l r I 1 -21 

+ o{n 'J) 



^2.+ir(s) 



h 

error2 = lim —ai(IC)By(^) 



Therefore, 



£^ = E[ijjK] + AMk + o(n') ^ £;[V'k] 



2v^ 



-| "0-1 



7[n — 1, v^] 



-,2n-2 



n=2 



+ 



H p^ai(/C) lim 



rrid 



2s+l 



7[s + i t)'] 



Note that the contributions proportional to ai(/C) of the poles at s = — 2 in Ai£-^(s) and A2£-^(s) 
cancel. 
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We are left with the very compact formula: 



AMk = hmd [Ao + Dn,] < 



Ao 
D 



no 



'2v^ 

"0-1 

-E 



n=2 



(9) 



57r t> 



2n-2 



In sum, there are only two contributions to semi-classical kink masses obtained by means of the 
asymptotic method: 1) /im^Ao is due to the subtraction of the translational mode; 2) hrridDno 
comes from the partial sum of the asymptotic series up to the no — 1 order. We stress that the 
merit of the asymptotic method lies in the fact that there is no need to solve the spectral problem 
of /C: all the information is encoded in the potential V{x). 



3 Loop kinks 

The existence of kinks is guaranteed if the minima of (/{tp) are a discrete set which is the union 
of orbits of the discrete symmetry group of the system. We shall use the term "loop" kinks to 
refer to those classical solutions that interpolate between vacua belonging to the same orbit of the 
symmetry group; otherwise, the solitary waves will be referred to as "link" kinks, see pnf. In this 
Section we shall discuss three kinks of the "loop" type. 

3.1 The quantum sine-Gordon soliton 

We first treat the sine-Gordon model by considering the potential energy density: U[ip{y'^)] = 
^(1 — cos ^ip). The dimensions of the parameters m and A are respectively: [m] = L^^ and 

[A] = M^^L^^. Therefore, we choose = m and Q = — and find: U[(f){x,t)] = (1 — cos0). 

The "internal" symmetry group of the system is the infinite dihedral group Dqo = Z2 x Z 
generated by internal reflections, — —(f), and 27r translations, — >■ + 27r. The vacuum classical 
configurations (l)\{x,t) = 27m form the orbit A4 = and there is spontaneous symmetry 
breakdown of the internal translational symmetry through the choice of vacuum. The moduli 
space of vacua, however, M. = is a single point and all the equivalents kinks of the model, 

4m 27mm 

0K(a;, t) = ±4 arctane + 27rn , -^kI?/, l/o) = arctane , 

V A V A 



are loop kinks. It is easy to check that -E['0k] = ^y- and -E[V'v] = 0. 
The second order variation operator around the kink solutions is 

(f - cP 

/C = — — + m^(l — 2 sech^ my) , K, = — — r + 1 — 2 sech^ x . 
dy^ dx^ 

Note that /C = im?K,\ henceforth, C,pk:{s) = ^^Cp/cI-^)- Simili modo, in the vacuum sector we have: 

V = -;^+m2 , V=-^ + l , V = m'V , Cv{s) = \M') ■ 
dy^ ax^ m^'^ 
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3.1.1 Exact computation of the mass and the wave functional 

• Generalized zeta function of V: 

The spectrum of V acting on functions belonging to L^(R) is Spec V = /c^ + 1, with A; G M a 
real number. There is a half-bound state fk'2=o{x) = constant that we shall not consider because 
it is paired with the other half-bound state in Spec(^). The spectral density on the interval 
/ = [— ^] with periodic boundary conditions is Pv(/c) = The heat function is, 



2vr y_oo 

and the generalized zeta function reads: 



Therefore, Cvi.^) (hence Cv('5)) is a meromorphic function of s with poles at s = |, — |, — |, — |, ■ ■ ■. 
The generalized zeta function of the Hessian at the vacuum is, however, also infrared- divergent: 
it is linearly divergent when L oo even at points s G C where (vi^) is regular. 
• Generalized zeta function of K,: 

In this case Spec^ = {0} U {k"^ + 1}, /c G M and the spectral density on I is 

mL 1 d6(k) 

with phase shifts 

6(k) = 2 arctan — 
k 

because fC is the Schrodinger operator that governs the scattering through the first of the "trans- 
parent" Posch- Teller potentials, ||T^. Thus, 

271- y_oo 27r dk 

= 1 + ^Le-'^ - Erfc 



where Erfc is the complementary error function, |18]. Note that /C has a zero mode, the 



eigen-function being the translational mode = sech x, which must be subtracted because it 
arises in connection with the breaking of the translational symmetry, x ^ x + a, by the kink 
solution and does not contribute to the kink mass up to this order in the loop expansion. There 
is also a half-bound state, fi.2^Q[x) = tanhx, that exactly cancels the contribution of the constant 
half-bound state in SpecV. Therefore, we obtain 

Tre"^-^^ = Tre"^^ - 1 = Tre~^^ - Erfc ^/p 

and 



oo 1 Y(s + 
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Cpk:("S) (hence Cpic{s)) is also a meromorphic function of s that shares all the poles with Cv{^)y but 
the residues are different except at s = a pole where the residues of (pici^) Cvi^) coincide. 
The infrared divergence, however, is identical in the kink background and the vacuum. 

We can now compute the limit of the regularized quantities that enter in the one-loop correction 
formula to the kink mass: 



K 



lim 



h 

2 
hm 

lim 

271 e^O 



711? 



hm 



limf^V' 



Tie) 



o\mJ r(i + e) 



i_21og^-V^(l)+V^(l) + o(e) 
em 2 



(10) 



and 
A,e^ 



2h 
T 

hm 



— — lim — 

L s^-l ^m 



lim ( — ) 



/i\2^+i r(s + 1) 



r(.) 
r(^) 



C^{s + i) + o{h^-f) = 

1 2/i 



hm 



Mm I iiV-" + 



TT s^-i \m/ 



T(s] 



vr e^o \mJ r(- 
r'(^) 



hm 

lim 

2n e~*o 



21og=^ + V^(l)-^(--; 
m 2 



oie 



+ o(n%) 



where ipiz) = is the digamma function. 

The important point to notice is that the renormalization of the zero-point energy performed 
by the subtraction of Cv(~|) still leaves a divergence coming from the s = — | poles because the 
residues are different. The correction due to the mass renormalization counter-term also has a 
pole. The sum of the contributions of the two poles leaves a finite remainder and we end with the 
finite answer: 



hm 



TT 



hm 



TT 



+ o{h^^) 



5m 

7 



hm 



TT 



+ o{h^^). 



'Ill 



The one-loop quantum correction to the mass of the sine-Gordon soliton obtained by means of 
the generalized zeta function procedure exactly agrees with the accepted result, see [Q, and, 
henceforth, with the outcome of the mode number regularization method, 0. 

The square of the modulus of the ground state wave functional up to one-loop order is given 
by (I). If = § and C = ^, obviously, Cw{s) = C^Xpds) (Ch^(O) = Cpyc(O)) and we have 

'^%). Thus, 



— - C^Xpk{s) log C + C^"^ 



ds 



PK 



ds 
ds 



dCv 1 r(g + |; 

ds y/TT sT{s) 

mL T(s — 



V'(s + i) 4j{s) 



ATI Vis) 



[^(3-|)-^(s)] 



;i2) 



and 

Cv(0) 
Therefore, 



Cpyc(O) = -1 



ds 



(0) = -mL , ^(0) = -mL + ^(l)-^(i). 



l*o(0K(x))r 



— exp i ^^-^ 



i*o(0v(x))r 



exp '^^L 
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Renormalizing the wave functional with respect to the vacuum we obtain 



3.1.2 The asymptotic expansion and quantum corrections 



(13) 



In the sine-Gordon model the exact formulas for Tre~^^'^ and (pjcis) are readily derived because 
the spectrum of the Schrodinger operator IC is completely known. On the other hand, the series 
expansion of the complementary error function tells us that 



1-3-5-...- 



l-3-5-...-(2n-l) 

and the a„(^) coefficients can be computed from this exact expression: 

, ao{IC) = mL , cin{fC) 



Tre 



-13 PK 



n=0 



i2n - 1) 



One can check by direct calculation that indeed, 

o-ni^-) = / dxan{x,x) ] n = 0, 1, 2, 3, . . . 



mL 
2 



and the a„(/C) are the integrals of the functions defined in Appendix for V{x) = 2sech^x. 

In any case we see from the formula (^) that the comparison with the exact result is satisfactory: 



710-1 



AMk = -0.282095^m - ^ a„(/C)rf„nm versus AMk = -0.318309;im . 

n=2 

The partial sums 

no-l no-1 r ^ 



n=2 



n=2 



can be estimated with the help of the following Table, 



n 


a„(/C) 


no-1 




2 


2.66667 


2 


-0.0670702 


3 


1.06667 


3 


-0.0782849 


4 


0.30476 


4 


-0.0802324 


5 


0.06772 


5 


-0.0805373 


6 


0.012324 


6 


-0.0805803 


7 


0.0018944 


7 


-0.0805857 


8 


0.00025258 


8 


-0.0805863 


9 


0.00002972 


9 


-0.0805863 
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For instance, choosing no = 10 we find that Dio = —0.080586 and the correction obtained by 
means of the asymptotic expansion is: 

E[iIjk] + ^Mk = ^ - 0.36268inm + o{h^\) 
A 

In fact 

hm,^ ,1, „ / iM , hm .^/lErfc-v/^ \ 



2 



is almost the total error: 0.044372/im. The difference is: 

hm 



&io,/c(-i)~10-'^m. 



4^/^ 

Note in the Table that a„(K^) rapidly decreases with increasing n. 
3.2 The quantum \{4)% kink 

We now consider the other prototype of solitary waves in relativistic (l+l)-dimensional field 

theory: the kink of the A(0^)2 model. The potential energy density is: U[ip{y^)] — | ~ x') " 

We choose, however, rud = ^ but keep Cd= ^ and find: U[(j){x,t)] = | (0^ — 1)^. 

The internal symmetry group is now the Z2 group generated by the — > — refiections and 
the orbit of vacuum classical configurations 0y(,x, t) — ±1 e M gives rise to a moduli space of 
vacua M = 1^ which is a single point. The kink solitary waves are thus loop kinks and read 

0K(a;,t) = itanha; , ^k(|/, 2/o) = tanh — . 

The kink and vacuum solutions have classical energies oiE[il)K\ = and E[iI)y] = respectively. 
The Hessian operators for the vacuum and kink solutions are 

^ ^ "d^^ ^ "^^^h^^T^rf^^ ^^^J ^ 2 
and the corresponding generalized zeta functions satisfy 

Cpk:(«)=|— ) Cpk{s) , Cv(«)=f— 1 Cvis). 
ml \ m 
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3.2.1 Exact computation of the semi-classical mass and wave functional 

• Generalized zeta function of V = — ^ + 4. 

Acting on the L^(]R) ® C Hilbert space we have that Spec V = {fc^ + 4}, A; G M, but the spectral 
density on the interval I = [—^^, of eigen-functions with periodic boundary conditions is 
Pv(^) = 273^" these data, the heat and generalized zeta functions are easily computed: 



Tre-^^ 



mL 



2V2' 



(vis) 



mL 1 

v^r(s) 



mL 



-4f3 



mL 1 r(s - I) 



(14) 



v/87f4"-5 r(s) 

and we find that C-]){s) has the same poles and infrared behaviour in the X{4>^)2 and the sine-Gordon 
models. 



Generalized zeta function of K, 



+ 4 V 

dx'^ cosh X ' 



JC is the Schrodinger operator for the second transparent Posch- Teller potential, ^7^. Thus, 
Spec IC = {0} U {3} U {k'^ + 4}, A; G M, and the spectral density on I is 



mL 



1 d5{k) 



where the phase shifts are 6{k) 
Tre-^^^ 



-2 arctan 



2V27r 2tt dk 

if PBC are considered. Thus, we find 



mL 



Stt 



O 1 

27r 



^dk ^ffie-'^(^'+^) 



mL 



+ e-=^^(l 



Erfcv/^) -Erfc2v^. 



The Mellin transform immediately provides the generalized zeta function: 



Cpds) = Cv{s) + 



T2Fi[\, S + i, |, -\] 



1 1 

4^ s 



(15) 



where 2-P'i[a, 6, c, d] is the Gauss hypergeometric function, |]T8 
The power expansion of 2-^1, 



2' 2' 



f. (-i)'r(/ + |)r(. + / + |: 



r(|) 



r(|)r(^+2M^o 



r(/ + |; 



tells us that, besides the poles of Cvl-s), Cpici^) has poles at s 



/ G Z+ U {0}; i.e., as in the sG soliton case, Cvl-s) and C,pk.{s) share the same poles except s = 2 
but the residues in the A(0^)2 model are increasingly different with larger and larger values of 
iRe s|. 
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Applying these results to the kink mass formula, we obtain 
Ai£^ -- 



2 \m'^ 



hm , f2fx^Y Tie] 
lim 



L2' 2' 3J 



— ■;=— lim 
2v/27r 



--31n^ + 2 + ln--2F^[|,0,|,-i] + o(5) 



hm — hm -^TT^^vis + 1 

L^oo L s^-i J T[s\ 



3hm, f2u'^Y A-'T(e) 
hm ' ' 



27re-*oVm2/ r(-i+£) 



3hm 

— ^ hm 

2V2tt 



2/i 



1, 



+ ln^ -ln4 + ?/;(l) - + o(e) 



e ' ' ' 2' 

where ^F'x is the derivative of the Gauss hypergeometric function with respect to the second 
argument. Therefore, Aie^ + A2e^ = ^ - ^ , and we obtain: 

4. - E\M + AMk ^ 1^ . ( - . o(S^A, , 

the same answer as offered by the mode- number regularization method [0. 

To compute the norm of the ground state wave functionals we closely follow the procedure 
applied in sub-Section §3.1. to the sine-Gordon soliton. In the A(0^)2 model, we find that 



PK 



dCv ^ 1 r(5-F|) 



4-^(- + ln4 + V^(s)-V^(s + |) ) - 



ds ds a/tt r(s) 

-2s 3-^-^ 2i^i[i s + i, |, -|](log3 - ^{s + \)+ ^{s)) + 2s 3"^^^2i^i'[i s + 



1 3 _ll 

2' 2' aJ 



and 



(^(,_ i)_^(,)_log4) . 



ds V87r4"-5 r(s^ 
from these expressions and formulas ([1^) and (|15|) one checks that 



(is 



Cv(0) = , Cp£(0) = -1 . 
We obtain 

The quotient of the probability densities is 

|*o(0K(a;) 



^(0) = -v^mL, ^^(0) = -y2mL + log48 





mL 


j exp 


.2^2. 



(is 



|^o(0v(a;))|^ = exp 



niL 
2^2. 



1 / C \ 2 



|*o(0v(x))|' 2 Vv^ 



(16) 
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3.2.2 The asymptotic expansion and quantum corrections 

In the X{(f)'^)2 model ^^\(i,y{x) = 4 and V{x) = ^§-|<v,v(x) — ^^Icpj^ix) = 6 sech^ x are the potentials 
of the Schodinger operators that respectively correspond to the Hessians at the vacuum and the 
kink configurations. The asymptotic expansion of the heat function 



-4/3 °° j-f^ -4(3 «^ 

Tre-^^'^ = -! + -== V / a,(a;, = -1 + V a,(£) 



n=0 " 2^2 

can be either obtained as a series expansion of the exact result 



IVe-'^^^ = -1 + 



(2n-l)!! 

or from the coefficients defined in the Appendix for V{x) = Gsech^a; 



-4p 



a„(/C) 



rriL 

dx 

mL 
2V2 



ao(/C) = , a„(/C) = 



V2 



(2n- 1)!! 



To compare with the exact result, we apply the formula given in the Appendix and observe 



that 



4 m' 



no-l 



^ ^ --^ - 0.199mhm - V an{IC)dnhm + o{h^\) 



is far from the exact result 



4 



3^A 



0.471imm + o(/i^A) 



before adding the contribution of the terms between n = 2 and n = no — 1 in the asymptotic 
expansion to the contribution coming from the subtraction of the translational mode. The partial 
sums 



no-l no-l r _ ^1 



n=2 n=2 

can be estimated up to no = 11 with the help of the following Table 



n 


a„(/C) 


no-l 


Dno 


2 


24.0000 


2 


-0.165717 


3 


35.2000 


3 


-0.221946 


4 


39.3143 


4 


-0.248281 


5 


34.7429 


5 


-0.261260 


6 


25.2306 


6 


-0.267436 


7 


15.5208 


7 


-0.270186 


8 


8.27702 


8 


-0.271317 


9 


3.89498 


9 


-0.271748 


10 


1.63998 


10 


-0.271900 
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Choosing no = 11, we find that Dn = — 0.271900/?-m and the correction obtained by adding 
Diihm is: 

AMk = -0.47137inm + o{h^X) 
in good agreement with the exact resuh above. In fact 



^[5p.(-|)-i?v(-|)] + ^5v(|) 

= —= / df3 ^ + + 3^ + ^ 0.00032792nm 

2v^ Ji V 2/32 2(3-2 2(32 ; 

is almost the total error: 0.0002580/lm. The deviation is 



hm , , 1 



4V27r 



With respect to the sine-Gordon model there are two differences: a) in the A(0^)2 model the 
error committed by using asymptotic methods is smaller, of the order of 10~^/im, a 0.07 percent, 
as compared with 10~^^m, a 6.00 percent, in the sG case; b) the rejection of the contributions of 
the no > 11 terms and the non-exact computation of the mass counter-term contribution has a 
cost of approximately 10~^^m in the A(0^)2 model versus 10~^^m in the sG system. Both facts 
have to do with the larger value of the smaller eigenvalue of the vacuum Hessian in the X{(f)^)2 
model with respect to the sG system, 4 versus 1. 

3.3 The cubic sinh-Gordon kink 

We shall now study a system of the same type where the potential energy density is: U[tlj{y'^)] = 
^ ^sinh^^'?/' — ij . Non-dimensional quantities are defined through the choice = m and 
Cd = the Euler-Lagrange equation is 

1 



n(f){t,x) = --sinh(20)(sinh> - 1) (17) 

and the justification for the choice of name is clear. We find this model interesting because it 
reduces to the X{4>'^)2 system if \(f){t, x)| < 1 and is the Liouville model, |T^, with opposite sign of 



the coupling constant, in the 4>{t,x) = ±cx) ranges. In fact, the potential energy density = 
|(sinh^0 — 1)^, see Figure 1(a), presents two minima at the classical values: 0v = ±arcsinhl. 
The two vacuum points are identified by the — > —(f) internal symmetry transformation and the 
semi-classical vacuum moduli space is a point. For this reason, U{x) has been apphed to the study 
of the quantum theory of diatomic molecules: the solutions of the associated time-independent 
Schrodinger equation are a good approximation to the eigen-states of a quantum particle that 
moves under the influence of two centers of force. We deal with the C = 1 ^ind M = 3 member of 
the Razavy family of quasi-exact ly- solvable Schrodinger operators, |2^, although we are looking 
at it from a field-theoretical perspective. 
The solutions of the first-order equations 

d(j) 1 / 2 \ tanh(x + b) 

— = ±—=(smh. (p — 1) , (I)k(x) = ± arctanh ■= , (18) 

dx y/2 V2 
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Figure 1: Graphic representation of (a) the potential energy density (h) the kink and (c) the Hessian 
potential well. 



see Figure 1(b) for 6 = 0, are the kink solitary waves of the system. The Hessian operators at the 
vacuum and kink solutions are respectively: 



/C 



V 

16m^ 



+ Am 



14m^ 



m 



m 



+ 4 



+2+ 



16 
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dy^ ' ' (1 + sech^ l + sech^y dx"^ (l + sech^x)^ l + sech^x' 

The mass of the fundamental mesons is thus 2m. K, is an Schrodinger operator: 



m K, 



rf2 



dx"^ 



+ 4 - V{x) 



V{x) 



2 sech^x(9 + sech^x) 



(1 + sech xY 

where the potential well plotted in Figure 1(c), albeit analytically very different from the s-G and 
A(0^)2 kink potential wells, exhibits a similar shape. 

We shall not attempt to solve the spectral problem of K,. The only thing that we need to know 
in order to apply the asymptotic method is that the lowest eigen-state is the unique zero mode: 

2V2 



dx (3 + cosh 2x) 

Therefore, the energy of the semi-classical kink state is approximately (see formula (^) 



m^ ( 1 "0^^ 

r\ ^ ^ (^1 _ 3V2arcsinhl) - Km — = + ^ a„(£)- 

\ * n=2 



87r4 



n-l 



(19) 



In the Table below we write the Seeley's coefficients and the partial sums = — Yln°=2^ Qra(A^) §^"4^-1^ 
up to no = 11: 



n 


a„(/C) 


no - 1 




2 


29.1604 


2 


-0.20135 


3 


39.8523 


3 


-0.26501 


4 


42.1618 


4 


-0.293253 


5 


36.0361 


5 


-0.306715 


6 


25.7003 


6 


-0.313005 


7 


15.6633 


7 


-0.315779 


8 


8.3143 


8 


-0.316917 


9 


3.9033 


9 


-0.317349 


10 


1.6590 


10 


-0.317502 
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obtaining the approximate answer: 

AMk = hm[Ao + Dii] = -0.282095nm - 0.317502nm = -0.599597nm 

We cannot estimate the error but we assume that this result is as good as the answer obtained 
for the A (0^)2 kink because the continuous spectrum of fC also starts at 4. 



4 Link kinks: the A((/)^)2 model 

Finally, we consider the following potential energy density: U[ip{y^)] = 4^'?/'^ (^i^^ — ■ The 
choice of rud = ^ and Cd = ^ leads to the non-dimensional potential: U[(f){x, t)] = (0^ — 1)^. 
The moduli space of vacua M = made out of two Z2 orbits, contains two points: 

Quantization around the (pvoi^: t) vacuum preserves the —>■ —0 symmetry, which is spontaneously 
broken at the degenerate vacua 0y^(a;, t). The kink solitary waves of the system 



1 / TTl / Tfl ~ 

0K(x,t) = ±-^Vl±tanh(a; + &) , ^K(y, = ±^ W 1 ± tanh-^(?/ + 6) 

interpolate between 0yj^(x,t) and 0y(,(a;,t), or vice- versa, which are vacua belonging to distinct 
Z2 orbits: these solutions are thus link kinks. 

The kink and vacuum solutions have classical energies of -E['0k] = -E[^/'yo] — EH'v±] = 

respectively. The Hessian operators for the vacuum and kink solutions are 

V, = -4^ + 2m^ = ^^44 + 4U^V. 



(P bm? 3m^ my Ibm? 
JC = — —7 H ± — ; — tanh^^ — 



dy^ A 4 V2 Scosh'^ 

mV 5 3 , 15 \ m^- 

= — — — H — ± -tanhx = — 

2 \ dx^ 2 2 4cosh2xy 2 

The problem of the semi-classical quantization of these and other link kinks have been ad- 



dressed somewhat unsuccessfully in ||21[] due to the analytical complexity of the eigen-functions of 
K, as well as the conceptual difficulty of dealing with a QFT on the real line where the asymptotic 
states far on the left and far on the right correspond to mesons with different masses. This issue 



has been analyzed in depth in |]22[: the main suggestion is that the normal-order prescription 
should be performed with an arbitrary mass to be fixed in order to avoid the ambiguity induced 
by the step function background. We now apply the asymptotic expansion of the heat function 
method in this complex circumstance to find a very natural way of choosing the mass renormal- 
ization parameter. Moreover, we improve the approximation obtained in the computation of the 
quantum kink mass by going farther than first-order in the asymptotic expansion. 
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Figure 2: Graphical representation of (a) the potential energy density (b) the kink and (c) the Hessian 
potential well. 



Besides the bound state, 



2cosh xa/2(1 ± tanhx) 



the u;^ = translational mode, the spectrum of K, includes transmissionless scattering states for 
1 < cj^ < 4, and states with both non null transmission and reflection coefficients if uj"^ > 4. In 
the language of QFT, the topological sectors based on link kinks are peculiar in the sense that 
the iV-particle asymptotic states are mesons that have different masses at a; = ±cx). If the meson 
energy is less than 2m^, the bosons are reflected when coming from the left /right towards the 
kink. More energetic mesons can either be reflected by or pass through the kink. If the mesons 
are transmitted there is a conversion from kinetic to "inertial" energy, or vice- versa, in such a way 
that the poles of the propagators far to the left or far to the right of the kink can only occur at 
= ^ and = 2w?. 

This is the reason why the subtraction from the Casimir energy of ^k, ^Cpk.{.~\) i of either the 
Casimir energy of the 0^^, |Cvo(~i)5 fCv±(~^)) vacua is hopeless, even after adding the 

mass renormalization counter-term to the Lagrangian. Therefore, we cannot use the generalized 
zeta functions Cvol-s) and Cvil-s) to renormalize the zero point energy in the kink sector. Instead, 
we will gauge the kink Casimir energy against the Casimir energies of a family of background field 
configurations that satisfy: 

501 (x) - 40|(a;) = ^(1 ± tanhax), (20) 
where a G M"*". The rationale behind this choice is that the a — oo limit is the background used 



by Lohe, ||21[|: 0^^(a;) = ±^(a:). The problem with Lohe's choice is that the discontinuity at the 
origin poses many problems for the algorithm of the asymptotic expansion because a nightmare 
of delta functions and their derivatives appears at a; = at orders higher than the first. Thus, we 
need some regularization, which is achieved by replacing the sign function by tanh in the formula 
(pOD above. In Figures 3(a) and 3(b) the Hessian potential wells for the backgrounds 0goo ^^'^ 



are compared. 

For any non- zero finite a, 05^(x) interpolates smoothly between ±| and ±1 when x varies 
from — oo to oo. The jump from ±1 to occurring at a; = in 0^^(a;) becomes a jump from 
±1 to 0, which therefore takes place at a; = ±oo!, followed by the smooth interpolation to ±1 

. If a = the background configuration is also pathological: 0ijo(a;) = ± — ^— ^,Va;, except at 
X = ±oo, where there are jumps to and ±1. 
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Figure 3: Graphic representation of the potential well produced by (a) the (f)Bi background (b) the (f)B, 
background as functions of x 



The Schrodinger operators 

» 5 , 3^ , 

"a — — rrr; H — ± -tannax 
dx^ 2 2 

« ^'5,3,, ^ 5 

^- = -^+2^2^^^) ' ^° = -d^ + 2 

govern the small fluctuations around the background (p^^. Thus, 

2 \ «+l 



A,s^(a) = ^ hm^ (^)^ ' [CpH^) - Cb„(«)] 



is the Casimir kink energy renormalized with respect to the (p^a background. From the asymptotic 
expansion of both Cpic{s) and CbA^) ^® obtain: 

where c„(^a) = ^ni^) " cin{Ba)- The deviation from the exact result is: 



hrud 

errori = 



In order to implement the mass renormalization prescription, we assume that virtual mesons 
running on the loop of the tadpole graph have a mass of ^ half of the time and a mass of ^/2m 
the other half-time on average. The normal-order is thus prescribed for annihilation and creation 
operators of mesons with M — mass; this amounts to considering 

as the infinite quantity associated with the single divergent graph of the system. Zeta function 
regularization plus the asymptotic expansion tell us that the induced counter-term adds 



=Ci(/Ca) lim 



"\-.-^,\5mlJ r{s) 
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to the one-loop correction to the hnk kink mass, whereas the error is 

errors = hm -^Ci(£«)5g„(|) . 
The sum of the contributions coming from the s — — | poles of Aie^{a) and A2e^{a) vanishes: 



lim 



7[s + 



- -1 

2' 2i 



r(s) 



7[g + i|] 

T(s) 



The choice of M = as a mass renormalization parameter leads to exactly the same result that 
we encountered in the more conventional systems with loop kinks and we end with the answer: 



AMk^ = hm[Ao + Dn,{a)] 



where An 



1 



2V27r 



and 



riQ-l 



no-1 



^no(a) 



Cn{lCa)dn = " ^ C„(/Cc 



n=2 



n=2 



n-1 



l[n - 1, 



The coefficients and the partial sums up to no = 11 for a 



Sv^TT 

1 are shown in the following Table 



n 


Cn(/Ci) 


no-1 


^no(l) 


2 


-9.3750 


2 


0.0968454 


3 


10.9375 


3 


0.0617547 


4 


-10.2567 


4 


0.0786049 


5 


7.89397 


5 


0.0703349 


6 


-5.12392 


6 


0.0741904 


7 


2.86874 


7 


0.0725233 


8 


-1.40987 


8 


0.0731872 


9 


0.61636 


9 


0.0729439 


10 


-0.24186 


10 


0.0730259 



We find: 



AMk^ ^ /im[Ao + Dii{l)] = -0.19947i;im + 0.0730259;im = -0.126445;im 



as the approximation to the kink Casimir energy measured with respect to the Casimir energy of the 
(j)§_^{x) background field configuration. 

The choice of a = 1 is optimum in the sense that for smaller values of a a tendency of the quantum 
correction towards — oo is observed whereas for a greater than 1 the tendency is toward -|-oo. In Figure 
4 , a = 1 is identified as the inflexion point of a family that interpolates between two background 
configurations with bad features: too abrupt if a = oo and too smooth if a = 0. 

We end this Section by comparing our renormalization criterion with the prescription used in [ p^ . 
Lohe and O'Brien choose a mass renormalization parameter M' in such a way that the mass counter-term 
exactly cancels the difference in vacuum Casimir energies between different points in the vacuum moduli 
space. 



V2 







(21) 
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a 



Figure 4: Quantum correction to the kink mass as a function of a in the (0.4, 2.0) interval 
The contribution of the tadpole graph must be considered for mesons with a suitable mass to satisfy (pl|): 

and we find M'^ = 2.33, M' = M';^, a very close value to M. At the L ^ oo limit 

CB-oa)-Cv(^) = iFiog|i 

If we had used M' as the mass renormalization parameter, the result would differ by 

,w A ,^ ci(^i) , 2.33 

Ak,M' - Ak,M = ^ ^ ;/ log ■ 



4^/2 (4 



TT 2 



2.50 



which is a very small quantity indeed. 



5 Outlook 

The natural continuation of this work, and the main motivation to develop the asymptotic method, 
is the computation of quantum kink masses in theories with N-component scalar fields. Nevertheless, 
explorations in the supersymmetric world along these lines are also interesting. 

All the models that we have described admit a supersymmetric extension because the potential energy 
density always can be written as U{tjj) = In non-dimensional variables the superpotential W{(p) 



for each model is: 



W{(l)) = ±4cos| 



W((/>) = ±4-^(isinh20-| 



The supersymmetric extension includes also a Majorana spinor field: 



a = 1,2. 



Choosing the Majorana representation 7^ = o" , 7 = ia = a of the Clifford algebra {7'^,7'^} = 2g' 
and defining the Majorana adjoint x = X^l^^ the action of the supersymmetric model is: 



S 



1 



dx^ { + ixi^df^x 



dW dW 



_d'^W 
X iio X 
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The M = 1 supersymmetry transformation is generated on the space of classical configurations by the 
Hamiltonian spinor function 

Q = J dx!^j^^\d^cP + n\^^ . 
The components of the Majorana spinorial charge Q close the supersymmetry algebra 

{Qa, Qp} = 2(7'^7°)a/3^^ - 27^/3^. (22) 

Their (anti)-Poisson bracket is given in (^) in terms of the momentum and the topological central 
charge T = \f dW\ . 

The chiral projections Q± = Q and x± = X provide a very special combination of the 
supersymmetric charges: 

+ Q- = / <ix {(x^ - x-)f - (X- + . 

+ is zero for the classical configurations that satisfy ^ = and x± = which are thus classical 
BPS states. One immediately notices that our kinks are such BPS states and besides the small bosonic 
fluctuations one must take into account the small fermionic fluctuations around the kink for computing 
the quantum correction to the kink mass in the extended system. The fermionic fluctuations around the 
kink configuration lead to other solutions of the field equations if the Dirac equation 

is satisfied. We multiply this equation for the adjoint of the Dirac operator 

f d'^W 1 f d'^W 1 

[-h^d^ + ^(<^k)| |i7^5^ + ^('Ak)| SFx{x,t) = 

and, due to the time-independence of the kink background, look for solutions of the form: 5p{x,t) = 
a;)e*'^*. This is tantamount to solving the spectral problem 

r d^ , d^W.^ ,d^W^^ , . ^dW,^ ,dm.^ , \ „ , , 2, / ^ 
[d^ ~d4P^^^^'~d^^^^' ^ ~d4^^^'~d^^^^'] '^^^''^^ " ^ fF{x;uj) . 

Projecting onto the eigen-spinors of ^7^, 
we end with the spectral problem: 

r ^d-'w d^w dw d^w \ M) (1) 2,(1), , 

for the same Schrodinger operator as that governing the bosonic fiuctuations. 

Therefore, generalized zeta function methods can also be used in supersymmetric models for comput- 
ing the quantum corrections to the mass of BPS kinks. Great care however, is needed in choosing the 
boundary conditions on the fermionic fiuctuations without spoiling supersymmetry. We look forward to 
extend this research in this direction. 
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Appendix 

In this Appendix we describe the iterative procedure that gives the coefficients a„(x,rc) used in the 
text. For alternative descriptions, see |23], |24|. For an interesting interpretation of these coefficients as 
invariants of the Korteweg-de Vries equation, see [26|. 

Starting from formula (P) in the text, we write the recurrence relation 



(n + 1) ttn+iix, y) + {x- y) - V{x)an{x, y) - ^ 



dx 

In order to take the limit y ^ x properly, we introduce the notation 

= lim^^^ 



dx"^ 



(23) 



and, after differentiating ( p3D k times, we find 

1 



^^^An{x) 



n + k 



{k- 



j=0 ^■'^ 



(x) 



from this equation and ('^'^^0(3;) = lim^^^a, = d'^^, all the ^^'^An{x) can be generated recursively. 
Returning to (p3|), we finally obtain a well-defined recurrence relation 



an+iix,x) 



n + l 



suitable for our purposes. 

We give the explicit expressions of the first eight a„(x,x) coefficients. The abbreviated notation is 



Uk = ^{x), ul - , 



df'V I 



ai(x, X 
a2{x,x 
a3{x,x 
04(3;, X 
05(2;, re 



no 
1 



1 



1 o 1 1 2 1 

—Ur, H UoUn H Ul H Ua 

6 ° 6 12 ^ 60 



^^0 + Y^-2no + -u,uo + -u,uo + ^u, + -u,u; + —ue 



11 



Ur, H UoUn H UtUr, H UMr, H Untin H UTU-iUt) H UfiUt) H UtU2 

120 ° 36 ^ ° 24 ^ ° 120 ° 40 ^ " 30 " 840 " 360 ^ 



26 



23 2 19 1 1 

+ 5040"^ + 2520™ + 280^^"^ + 15120"« 



1 4 1 61 o 43 23 2 5 2 19 

+ "Ui H usun H Ur, H U1U2U3 H uoUn H UtUa H unuoUA 

288 ^ 15120 15120 ^ 2520 ^ ^ ^ 5040 1008 ^ 2520 

23 2 19 1 2 11 1 1 

+ 30240^^ + 15120^=^"^ + 1680^0"^ + 15120^^^^ + 3780^^^^ + 332640^^° 

^) = 5M0 "0 + T^^^^o + i^'"^ + 2io"^^o + io + ^^'^2^0 + 5^^^^^"o 

I 4 61 o 43 5 2 1 23 2 2 
+ 288^^"° + 15120"^"'^ + 2520"^"^"^"° + 1008^^"^"'^ + 33264o"^°"° + 1008o"^"° 

19 2 1 ^ 83 2 2 1 31 2 1 1 4 

5040 ° 5040 ° 10080 ^ ^ 252 ^ 10080 ^ 280 1440 ° 

5 2 23 2 1 19 71 2 1 2 

+ 2016^^"^ + 30240"°"^ + 420"^"^"^ + 15120"°"^"^ + 665280^^ + 2016"^"^ 

II 61 1 19 1 2 17 

+ UnUoUe H UaUq H Ur\UiU7 H 1*3^7 H UnUa -\ UoUg 

15120 332640 3780 166320 30240 ° 332640 

1 1 



66528 ^ 8648640 



1 6 17 4 83 2 2^ 61 2 3 , 1261 4^ 43 2 

4320^°^^ + 8640^^^^ + 10080^°"^^^ + 30240^°^^ + 1814400^^ + 5040^°^^^^^=^ 

227 2 23 o 2 659 2 2 5 2 2 19 s 527 2 

37800 ^ 2 d -r 3Q24Q 3 ^ 302400 ^ ^ 20I6 " ^ ^ ^ 15120 ^ * ^ 151200 ^ 

7939 2 6353 2 1 s 17 s 13 3067 

9979200"^"^ + 9979200"^"^ + 1680^°"^"^ + 30240"^"^ + 12320"^"^"^ + 498930o"^"^"^ 

71 2 1 4 3001 2 13 61 

unUr^ H UnUe H Ur,ue H uiUsUr H uqUaUg 

665280 ^ 20160 ° 9979200 ^ 29700 332640 

3433 2 109 19 1501 71 2 

259459200 ^« + 498960^^"^^^ + 166320^°"^^^ + 64864800^^^^ + l99584o"^^« 
17 2003 1 5 1 2 

UnUoUg H UAUg H UnUiUg H UrUo H ^^n^tlO 

332640 129729600 66528 648648 665280 ° 

73 1 1 

25945920^^^^° + 1441440^^"^^ + 259459200 ""^^ 



+ 
+ 

+ 
+ 
+ 

+ 
+ 
+ 
+ 
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